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Abstract. Let X be a projective scheme; let M and Af be two coherent Ox- 
modules. Given an integer m, we present an algorithm for computing the global 
extension module Ext m (X; A4, AT). In particular, this allows one to calculate 
the sheaf cohomology H m (X, Af) and to construct the sheaf corresponding to 
an element of the module Ext 1 (X; A4,Af). This algorithm can be implemented 
using only the computation of Grobner bases and syzygies, and it has been 
implemented in the computer algebra system Macaulay2. 



Introduction 

Let X be a projective scheme over the field k; let Ai and Af be two coherent 
Ox-modules. We write Ext m (X; Ai, Af) for the m-th global extension module. 
There are several equivalent definitions for these modules. Grothendieck |^) and 



Hartshorne [RH] define them as the right derived functors of Af *— ► Hom(X; Ai, Af), 
where Hom(X; .A/f, AT) denotes the group of Ox-module morphisms. Griffiths and 
Harris |G-HC define the m via hypercohomology and discuss the correspondence be- 



tween the equivalence classes of global extensions of Ai by Af, also called Yoneda 
Ext, and Ext m (X ; Ai , Af) . The global extension modules play a central role in 
duality theory — again see Hartshorne [RH| or Griffiths and Harris ]G-H ] for 



more details. The purpose of this paper is to provide an algorithm for comput- 
ing Ext m (X; M,Af). 

To achieve this goal, we must determine how to represent the coherent sheaves 
Ai and Af. Because we are working over a projective scheme, every coherent sheaf 
can be represented by a finitely generated graded i?-module, where R is the homo- 
geneous coordinate ring of X. By embedding X into P™, we may identify R with 
a quotient of the polynomial ring 5* = k[xo, . . . , x n ]. This representation is advan- 
tageous since it allows one to use the algorithms and techniques of computational 
algebra — see Vasconcelos jvl. However, there is no canonical choice for the finitely 
generated module associated to a given coherent sheaf. In particular, the natural 
candidate ij£Z T(X, A/"(v)) may fail to be finitely generated. With this in mind, 
our problem is more precisely stated as follows: given two finitely generated graded 
i?-modules M and N, calculate the i?-module Ext m (X; M,Af) where Ai and Af 
are the coherent sheaves associated to M and N respectively. 

The answer for m = motivates our general result. In this case, Ext (X; Ai, Af), 
which equals Hom(X; Ai, Af), can be identified with (Hom^(Af> r , -AT)) for r>0, 
that is the set of degree zero i?-linear maps from a truncation M> r to N, for 
sufficiently large r. For a simple example in which truncation arises, we consider 
M = R and N — m d with d > 0; the second module is the d-th power of the 
irrelevant ideal in R. Both of these modules correspond to the structure sheaf Ox ■ 
Since a morphism from Ox is determined by the image of the unit global section 
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and T(X, O x ) = K it follows that Hom(A; O x , O x ) — k. Analogously, when r > d 
we have (Hom^(i?> r , Tn d )) Q = k. However, for r < d there are no degree zero maps 
from R> r to m d . Hence, to evaluate Ext (A; A4, AT) using this bijection, we must 
determine an appropriate truncation of M. To provide an algorithm, we need an 
effective bound r such that Ext°(X;M,N) = (Hom fi (M> r , N)) for all r > r Q . 

Our solution to the general question involves relating Ext m (A; Ai, Af) to the 
R- module Ext^(M> r , N); the second Ext module is calculated in the category 
of graded i?-modules. This method has the added virtue that it allows one to 
understand ® veIj Ext m (X; A4,Af(v)) as an i?-module. To express our main result, 
we write s_/V for the S-module obtained from ./V by restriction of scalars, aj(<?iV) 
for the maximal degree of the minimal generators of the i-th syzygy module of s A, 
and pdg(giV) for the projective dimension of gA. 

Theorem 1. Let M and N be two finitely generated graded R-modules, let m be 
an integer and set I = min{dirn(A), to}. If r is an integer satisfying 

r > max {a,i(sN) — i: n — £<i< pd s (siV)} — to + 1, 

then we have an isomorphism of graded R-modules 



e 

v>0 



Ext m (X;M,M(v)) S (Ext^(M> r ,7V) 



In particular, when pd s (s^V) < n — t, we have max{0} = — oo and the inequality 
for r is vacuously satisfied. Furthermore observe that, by replacing N with -/V(e) 
for some integer e, one can compute the module ©„ >e Ext m (X, Ai, Af(v)). We also 
point out that the integer r is closely related to the Castelnuovo-Mumford regularity 
of s-/V. This theorem generalize Eisenbud's approximation method for computing 



sheaf cohomology [E2, theorem 8.3.2]. 

We prove the main result in two steps. First, we examine the spectral sequence 
relating a locally free resolution T t of A4 to the global extension modules, namely 

=H p (ff 9 (X>m 0x (f„^))) ^Ext p+ «(X;X,AA). 

Applying vanishing conditions for sheaf cohomology, we obtain conditions for r 
which insure that spectral sequence becomes E™'° = Ext m (A; Ai, AT). The prob- 
lem therefore reduces to computing the cochain complex Hom(A; J 7 ,, A"). In the 
second step, we assume T» corresponds to a free resolution F, of M> r and consider 
the natural homomorphisms from Homj; [F^ , N) to ©„ eZ Hom(A; Ti, M(v)) . The 
kernel and cokernel of these maps are both local cohomology modules. By using 
local duality, we place further restrictions on r and bound the maximal degree of 
nonzero elements appearing in the kernel and cokernel. Combining both parts, we 
obtain the required isomorphism. 

The complexity of the algorithms derived from theorem [l] depends on computing 
Grobner bases. Because of significant differences between worst case bounds for 
Grobner bases and the complexity of geometric examples, we omit an analysis of 
the computational complexity. Instead, we refer the reader to Bayer and Mumford 
|B-M I for a discussion of the conjectures and results in this area. 



Background material is provided in the first section. In the second section, 
we prove two propositions which determine where to truncate M. A proof of 
theorem [l] is presented in the third section, followed by algorithms for computing 
global extension modules and sheaf cohomology. Section 4 contains four sample 
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calculations: an example in which the bounds in theorem [l] are sharp, a comparison 
of two methods for computing the cohomology of a locally free sheaf, an illustration 
of Serre-Grothendieck duality and the construction of the sheaf associated to an 
element of Ext 1 (X;M,Af). 

1. Preliminaries 

We collect here a number of definitions, standard results and notations. Refer- 
ences for the unproved assertions about free resolutions can be found in Eisenbud 



|E1|, the results from algebraic geometry are in Hartshorne [FtH], local cohomology 
is presented in Brodmann and Sharp p3-S| and the treatment of hypercohomology 
follows Weibel @. 

Free Resolutions. Let M and N be finitely generated graded i?-modules. If 
e is an integer, M e denotes the e-th graded component of M. The submodule, 
M> e = (§) d>e Md, consisting of elements of degree greater than or equal to e, is 
called a truncation of M. For an integer v, we write M(v) for the v-th twist of M 
defined by the formula M(v) e = M e+V . When we write isomorphisms and exact 
sequences of graded modules, we will always arrange that the maps have degree zero. 
However, the set of all homogeneous maps from M to N of all degrees is denoted 
Honifl(M, N); it is a graded i?-module, graded by the degrees of the maps. Using 
graded free resolutions, this construction extends to a grading on the R- modules 
Ext£(M, N). 

The (unique) minimal graded free resolution of M will be denoted F, (M) and 
pd R (M) will be the length of F, (M), in other words, the projective dimension of 
M. In particular, there is an exact sequence 

— » F p (M) — » > Fi(M) — » F (M) — » 0, (FR) 

where p — pd R (M). Here each Fi(M) is a direct sum of twists of R: 

bi{M) 

Fi{M) = R(- ai>j (M)) 
i=i 

and hence the maps in (FR) are given by matrices of homogeneous forms. In this 
setting, minimality means that none of the entries in these matrices are nonzero 
constants. 

We define Oi(M) to be the maximal degree of the minimal generators of the 
i-th syzygy module of M. Similarly, aJJSt) is the minimal degree of the minimal 
generators of the i-th syzygy modules. Rephrasing, we have 

ai{M) = max{a itj {M) : 1 < j < b t (M)} 

Q.i(M) = min {a itj {M) : 1 <j< h(M)}. 

We consider both a, (M) and (M) as functions defined for all integers i by using 
the conventions max{0} = — oo and min{0} = oo. We emphasize that cii(sN) is 
determined by a free resolution of s_/V in the category of S'-modules. 

Algebraic Geometry. By definition, we have X = Proj(i?) and P™ = Proj(S'). 
Sheaves of Ox-modules will usually be denoted by characters in script font, such 
as Q. The Ox-module associated to an i?-module M will be written M.. We point 
out that the free resolution F,(M) of the i?-module M gives rise to a locally free 
resolution !F,(M) of the sheaf M. The projective embedding of X is given by the 
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morphism of schemes /: X — > P™. As / is a closed immersion, the direct image 
functor /* takes coherent Ox-modules to coherent Op^-modules. Moreover, /*A4 
may be identified with the sheaf obtained from Ai by extending by zero. 

The global section functor will be written T(X, — ) and the q-th right derived 
functor, called sheaf cohomology, will be denoted H q (X,—). The sheaf of local 
homomorphisms, also called sheaf Horn, is denote 7iomo x (S, Q') and £xt q (Q, — ) 
is the q-th right derived functor of Ttomo x {G , — )• 

Local Cohomology. The maximal ideal of S, which is generated by the elements 
xq, . . . , x n , is called n. The image of n under the canonical map gives the maximal 
ideal m of R. For each _R-module M, we define T m (M) = lJ fcgN (0 :m Tn fc ); the set 
of element of M which are annihilated by some power of m. Notice that T m (M) is 
a submodule of M. The q-th right derived functor of r m (— ) is denoted by Hm(-) 
and is referred to as the q-th local cohomology functor with respect to m. The local 
cohomology modules Hm(M) are naturally graded Artinian R- modules. 

Hypercohomology. The p-th right hyper-derived functor of T(X, — ), also called 
the p-th hypercohomology functor, is denoted H P (X, — ). Let £* be a cochain 
complex of Ox-modules. We write H q (£') for the q-th cohomology group of the 
complex; that is the kernel of the q-th differential modulo the image of the (q — l)-th 
differential. For a given cochain complex, there are two spectral sequences relating 
the sheaf cohomology functors to the hypercohomology functors: 

'E p 2 ' q = W(H q (X,£')) =£» W +q (X,£ m ) (SS.l) 

"E%' q = H*(X,H. 9 (£ m )) =2* H p+q {X,£'). (SS.2) 

Both spectral sequences converge when £ * consists of only a finite number of 
nonzero Ox-modules. 



2. Bounds for Approximating Global Extension Modules 

We begin by giving vanishing conditions for the m-th cohomology group of a co- 
herent sheaf. This lemma can be viewed as refinement of Serre's finiteness theorem 
|RH| , theorem III. 5. 2] or Cartan's theorem B []G-H , section 5.3]. 

Lemma 2.1. Let N be a finitely generated graded R-module and let m be a positive 
integer. If the inequality v > a n -m(sN) — n is satisfied, then the cohomology group 
H m (X,Af(v)) vanishes. 

As hypercohomology will be required in the following discussion, we provide a 
proof using this machinery. One could also prove this assertion by induction on the 
projective dimension of giV. 

Proof. Given the closed immersion /: X — » P", there exists an isomorphism of 
cohomology groups [RH, lemma III. 2. 10]: 

H m (X,M{v)) S H m (¥ n ,fM{v)). 



As f*Af(v) is the sheaf associated to the S'-module s-/V(v) [RH, proposition II. 5. 12], 
it suffices to prove the assertion for projective space. 

Without loss of generality, we may assume X = P™, that is R — S. Consider the 
cochain complex given by £~ % = (^(N)) (v). As f»(N) is a resolution of A/", the 
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second hypercohomology spectral sequence (SS.2) collapses and the first spectral 
sequence (SS.l) becomes 



E p 2 ' q = R p (H q (X,£')) =£*■ H p+q (X,Af{v)). 



We recall that, for projective n-space, one has the following equivalence [RH, the- 
orem III. 5.1]: 




(VC) 



H*(X,O x {v))=0 
Because F m (N) is a free resolution of N, we have 

bi(N) 

= O x (v-a id (N)) 

for < i < pd R (N); otherwise £~' = 0. It follows that the relation E% q ^ 
implies q equals or n and ~pd R (N) < p < 0. From this, we see that, for p + q > 0, 
either the target or source of the differential d^' q : 'E^ q — ► 'E p . +r ' q ~ 7 + is zero. We 
conclude that H m (X,Af(v)) = W n - ,l (H n (X,£ m )). By making use of the vanishing 
conditions (VC) for a second time, we observe that, when — pd H (iV) < m — n < 0, 
the inequality v > a„_ m (iV)— n implies that we havei?"(X,£:"- m ) =0. Otherwise 
we have £ n ~ m = and the vanishing is immediate. Finally, our convention implies 
that a n - m {N) — — oo for m — n < — pd R (N) or m — n > and the assertion 
follows. □ 

Remark 2.2. We may give a sharper version of the preceding lemma by placing 



conditions on the integer m. In fact, Grothendieck's vanishing theorem [RH, the- 
orem III. 2. 7] shows that the cohomology group H m (X,Af(v)) is zero when m is 
greater than the dimension of the support of Af(v), that is to > dim(AT) . 

Example 2.3. We recall that, for a complete intersection X in P" of dimension d, 



one has H m (X, Ox(v)) = for < m < d a nd a ll integers v RH , exercise III. 5. 5]. 
This illustrates that the converse of Lemma |2.l| is false. O 

Applying the above lemma, we provide bounds which guarantee that the first 
hypercohomology spectral sequences (SS.l) collapses. In these circumstances, we 
are able to compute Ext m (V; M, Af(v)) as the cohomology of a given cochain 
complex. More explicitly, we have: 

Proposition 2.4. Let M and N be two finitely generated graded R-modules, let 
to be an integer and set I = min{dim(iV), to}. For any integer v satisfying the 
inequalities 

v > max \a n - u {sN) — o m _ u (M) : 1 < u < 1} — n, and 
v > max {a n ^ u+1 ( s N) - a m _ u (M) : 2 < u < l} - n, 
we have the following isomorphism of graded R-modules: 

Ext m (X;M,Af(v)) = H m (Hom(V; F.{M), N{v)) 
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Proof. For m < 0, both the global extension module and the cohomology group 
vanish. Thus, we may confine our attention to m > 0. Consider the cochain 
complex defined by £ l = TLom o x (j : i(M),M{v)) . Because !F,(M) is a locally free 
resolution of Ai, we have [R.H, proposition III. 6. 5]: 

H*(£*) = w(nom 0x {F.{M),N{v))) = £xt q 0x (M,Af(v)) . 

Furthermore, for any two Ox -modules Q, Q 1 ', Theoreme 4.2.1 in ]G| provides the 
spectral sequence HP(X,£xt q 0x (g,g')) Ext p+q (X; Q, Q'). Hence, the second 
hypercohomology spectral sequence (SS.2) becomes 



„ E P,q = H p 



(x,£xt q 0x (MM{v))) ^Ext p+q (X;M,Xf(v)) = W +q {X, £ '), 



and consequently the first hypercohomology spectral sequence (SS.l) is 

'E% q = w(H q (x,7iom 0x (M M )^(v)))) =^ Ext p+q (X; M,M(v)). 
We point out that 

H (x,H)m Ox {F.{M),N{v))} =Bom(X;F.(M),Af(v)), 

whence the relation El? = W (liom(X; T.(M),Af(v)) 

The proof therefore reduces to showing that the bounds on v imply that the 
first hypercohomology spectral sequence (SS.l) degenerates; more precisely, we will 
prove 

p+q—m 



To achieve this, it suffices to demonstrate that 'E 2 ' vanishes for 1 < u < m and 
, E m-u,u-i vanishes for 2 <u <m, since the differential d p ' q : Ep q -> £»+n?-r+l 
has bidegree (r, 1 — r) and E p ' q — when q < or p < 0. Now, each !F p (M) is 



locally free of finite rank, so nom 0x (^ P (M), Af(v)) S 0^ ( X M) Af (v + a p ,j (M)) , and 
we obtain 

'b.(M) 



'El' q = W>\ iJ 9 (x,jV(« + a., J (M)) 



From remark 2.2, we immediately see that g > dim(./V) implies 'E 1 ^' 9 = 0. Moreover, 
applying lemma |2.1| , we see that if the inequality v > a n - q (sN) — a p (M) — n is 
satisfied then 'i?^' 9 is zero. Combining these observations, the claim follows. □ 

Next, we relate the m-th global extension module to the graded i?-module 
Ext£(M, N). 

Proposition 2.5. Let M and N be two finitely generated graded R-modules and 
let m be an integer. For any integer e satisfying the inequality 

e > m&x{a n+1 ( s N),a n ( s N)} - a m (M) - n, 

there is a natural isomorphism of graded R-modules 

(Ext£(M,iV)) >e £ 0H m (Hom(X;^.(M),AA( V )) 
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Proof. Again, for m < 0, both the global extension module and the cohomol- 
ogy group vanish, so we may assume that m > 0. Now, the graded i?-module 
Ext^ l (M, A) can be computed from a free resolution of M, that is 



Ext£ (M, TV) = H'» (Hom fl (F. (M) , 



N 



Relating local cohomology to the sheaf cohomology, we have the following natural 
exact sequence JB-S|, 20.4.4]: 



0^ H^RomRiFiiMlN)^ — > Hom fl (F i (M ), A) — ► 
-^0Hom(X;^(M),7V(w)) (Hom s (F i (M),A)) — .0, 

for each i; in other words, we have a chain map 

Homfi (F, (M) , A) — > 0Hom(X;jr.(M),Ar(u)). 

Thus, it suffices to show that this chain map is an isomorphism for degrees greater 
than or equal to e. 

Using the properties of the bifunctor Homfl(-, — ), we have 

(b m (M) \ fc m (M) 

R(-a mtj (M)),N\ = A(a mj -(M)). 
i=i / i=i 

Moreover, as ni? = m, the A:-th local cohomology module of the i?-module A, 
H^(N), i s iso morphic to the A;-th local cohomology module of the S- module 5 A, 
H^(sN) p^S| , 13.1.6]. It therefore suffices to prove that, for < k < 1, the 
maximal degree of an element in ffJ(giV) is less than e + a m (M). With this in 
mind, define Ck to be the maximal degree of an element in H^(sN). Notice that the 
local cohomology modules H^(sN) are Artinian, so ct is well defined. Applying 
the graded form of local duality |B-S|, 13.4.6], we have 



H%N) £* (Rom k (Ext^ +1 - k ( s N,S),k))(n + 1 



Setting c k to be the minimal degree of an element in Ext" +1 ~ fe ( s 7V, S), we obtain 
the formula Cfc = — c k — n — 1 . To compute a lower bound for c k , we use the free 

resolution F.( S N) of S A; in particular, Ext|( s A,5) = H fe (Hom s (F.( s A), 5)) . 

As above, the properties of Homs(-, — ) give 

bi( s N) 

Rom s (Fi( s N),S)^ SfajQiN)). 

Consequently, we have c k > —a n +i~k(sN) which yields the inequality 

max{co,ci} < max{a„(sA), a„+i(gA)} - n. 

The conditions on e clearly imply that e > max{co , ci } — a m (M) which completes 
the proof. □ 

We end this section by summarizing our results. 
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Corollary 2.6. Let M and N be two finitely generated graded R-modules, let m 
be an integer and set £ = min{dim(A r ), m}. If the integer e satisfies the following 
inequalities: 



77; 

n, 



e > max {a n ^ u ( s N) - a m _ u {M) : < u < £} 
e > max [a n - u+ i( s N) - a m _ u (M) : <u < £, u^l}- 

then we have an isomorphism of graded R-modules 

(Ext™ (M, iV)) Ext" 1 (X; M,M(v)) ■ 



Proof. This follows immediately from proposition p.4| and proposition 2.5.D 



3. A Method for Computing Global Extension Modules 

Using the bounds developed in the preceding section, we present a proof of the 
main theorem. 

Proof of Theorem Fix an integer r such that 

r > max {a,i(sN) — i : n — £ < i < pd s (s-ZV)} — m + 1. 
To simplify notation, set M' = M> r . We remind the reader that the Ojf-module 



associated to the R- module M is equal to the one associated to M' [REj, exer- 
cise II.5.9]. Applying corollary ^.6j it suffices to show that the following inequalities 
are satisfied: 

a m - u {M') + n-a n - u { s N) > for 0<m<^; 
a m _ u {M') + n-a n - u+ i{ s N) > for 0<u<£,u^l. {> 

By convention, a m _ u (M') = 00 for m-«<0orm-u> pd R (M), so we may 
confine our attention to m — pd R (M) < u < m. As Fq(M') maps onto M' and 
the maps Fi(Al') — > Fi^i(M') are given by matrices of homogeneous forms of 
positive degree, we have the inequality Oj(M') > r + i for all integers i. By making 
substitutions for u, we transform the inequalities (I) into the following: 

r > lli(sN) — i — m for n — £ < i < n; 
r > a,i(sN) — i — m + 1 for n — £+l<i<n+l, i =/= ji. 

Again, our conventions imply that ai(sN) = —00 for i < or i > pd s (sN), so we 
need only consider < i < pd s (siV). The choice of r therefore guarantees that 
these inequalities hold.D 



Remark 3.1. For dim(M) = 0, it is necessary and sufficient that M have finite 
length |El| , corollary 2.17]. In particular, if dim(Af) = then M has only finitely 
many nonzero graded components and we have M> r = for r > 0. It follows from 
theorem |l] that Ext m (X; M, N) = when dim(M)~= 0. 

For the special case X = P n , we have the following additional vanishing result: 

Corollary 3.2. Let M and AT be two coherent Ox -modules and assume X = P n . 
If m > n, then the global extension module Ext m (X ; A4 , Af) vanishes. 
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Proof. Let R be the homogeneous coordinate ring for X, and let M, N be two 
finitely generated graded R- modules which correspond to A4, J\f respectively. By 
theorem [j], we are required to show that, for m > n, the module (Ext^(M> r , AT)) 
vanishes for r sufficiently large. It therefore suffices to prove that pd R (M> r ) < n 
for r>0. Applying remark 3.1, we may assume dim(M) > 0. Now, since X = P™, 
we have R = S = k[xo, . . . , x n ] and the Hilbert syzygy theorem El, theorem 1.13] 
provides the inequality pd R (M> r ) < n + 1. Moreover, as R is a Cohen-Macaulay 
ring, we have ~pd R (M> r ) = dim(i?) = n + 1 if and only if the maximal ideal m of R 
is associated to M> r [El, corollary 19.10]. With this in mind, recall that H^(M) 
is the set of elements of M which are annihilated by some power of m. Since the 
module H^(M) is Artinian, there exists a maximal degree Co of an element in 
H^(M). Hence, if r > Co then the maximal ideal m is not associated to M> r which 
completes the argument. □ 

Using Grobner bases, one can compute the minimal free resolution of a finitely 
generated graded S- module (see chapter 15 in Eisenbud [El]). This allows one to 
calculate the projective dimension of a finitely generated graded S'-module and the 
degrees of the minimal generators of the syzygy modules. Furthermore, using the 
free resolutions one can construct the Ext modules for a pair of finitely generated 
grad ed i?-modules — again see chapter 15 in Eisenbud [El|. Grayson and Stillman 
G-S] have implemented these functions in their software system, Macaulay2. With 
this in mind, theorem [l] yields algorithm 3.3. By extracting the appropriate graded 
part, we have algorithm |3.4| which returns the m-th global extension module. For 
the special case M = R, we have the equality [HH, proposition III. 6. 4]: 

Ext m (X;O x ,M(v)) =H m (X,N(v)). (EH) 

Hence, we obtain the following two techniques for computing sheaf cohomology: 
algorithm 3.5 and algorithm 3.6. These algorithms have been implemented in 



Macaulay2; the routines may be found at ittp : //math.berkeley . edu/~ggsmith. 

The complexity of algorithm 3.3 will depend on computing F,(sN) and F,(M> r ); 
the second resolution is used to determine Ext^(M> r , N). To find these resolutions 
involves several Grobner bases calculations. As noted in the introduction, there 
is a significant difference between worst case bounds for Grobner bases and the 
complexity of geometric examples. With this in mind, we o mit an analysis of 
complexity and refer the reader to Bayer and Mumford | B-M for a discussion of 
the conjectures and results in this area. 



Remark 3.7. Algorithm |3.5| is equivalent to the appro ximatio n method for com- 
puting sheaf cohomology described in theorem C.3.1 in | E-G-S |. 



4. Some Examples 

All computations in this section were done on an Intel Pentium 166 MHz machine 
running Macaulay2 (version 0.8.41) for the Linux platform. 

We begin by giving an example in which the bounds in theorem [l] are sharp. 

Example 4.1. Let X = P 3 , let M = Ox, let J\f be the structure sheaf of the 
smooth ration al quar tic curve in P 3 and set m = 1. Following Eisenbud, Grayson 
and Stillman [E-G-S , these construction in Macaulay2 are: 

il : kk = ZZ/32003; S = kk[w,x,y,z]; 
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Algorithm 3.3. globalExtSum(m, e, M, N) 

input: two integers m, e and two finitely generated graded i?-modules M, N. 
assumption: the graded ring R is a quotient of a polynomial ring S. 
output: the graded il-module Q) v>e Ext m (X; M, Af(v )) , where X = Proj(_R) and 
Ai, N are the coherent Ox-modules associated to M, N respectively. 

begin 

if dim(M) = or m < then E := 
else ( 

n : = the number of generators of S minus 1; 
£ : = min{dim(./V), m}; 

$N := the S-modulc obtained from N by restriction of scalars; 

F,(gN) := the minimal free resolution of gN; 

pd( s N) := the length of F.( S N); 

if pd( s N) <n-£ then E := Ext™(M, N) 

else ( 

for i from n — £ to pd(gN) do ( 

cii(N) := the maximal degree of the minimal generators of Fi(gN); ); 
r := max {a;(sA r ) : n — i < i < pd(sAf)} — e — m + 1; 
B:=Ext™(M> r ,iV); ); ); 
e' := minimal degree of the generators of E; 
if e' < e then E := i?> e ; 
return E; 
end. 



Algorithm 3.4. globalExt(m, M, N) 

input: an integer m and two finitely generated graded R- modules M, N . 
assumption: the graded ring R is a quotient of a polynomial ring S. 
output: the k-vector space Ext m (X; M., M), where X = Proj(iJ) and M, Af are the 
coherent Ox -modules associated to M, N respectively. 

begin 

E := globalExtSum(m, 0, M, N); 
V := the k vector space isomorphic to Eq; 
return V; 
end. 



Algorithm 3.5. sheaf CohomologySum(m, e, N) 
input: two integers m, e and a finitely generated graded R-module N. 
assumption: the graded ring R is a quotient of a polynomial ring S. 
output: the graded R-module @ v>e H m (X,Af(v)j , where X = Proj(i?) and M is the 
coherent Ox -module associated to N. 

begin 

E := globalExtSum(m, e, R, TV); 
return E; 
end. 



Algorithm 3.6. sheaf Cohomology(m, N) 
input: an integer m and a finitely generated graded ij-modulc N. 
assumption: the graded ring R is a quotient of a polynomial ring S. 
output: the k-vector space H m ^X, Af(v)^J , where X = Proj(J?) and Af is the coherent 
Ox-module associated to N. 

begin 

E := globalExt(m, R, TV); 
return E; 
end. 
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13 : I = monomialCurve(S,{l,3,4}) 

3 2 2 2 3 2 

o3 = ideal (x*y-w*z, y -x*z , w*y -x z, x -w y) 

03 : Ideal of S 

14 : N = S'l/I 

04 = cokernel {0} I xy-wz y3-xz2 wy2-x2z x3-w2y I 

1 

04 : S - module, quotient of S 

Next, we determine max |a,(giV) — i: n — £<i< pdg(siV)}. Recall the minimal 
free resolution of N |E2| ; 8.1]: 

O^S(-5)^0S(-4)^S(-2)9 05(-3) — 5; 

i=i \i'=i / 

as n = 3 and I = 1, we obtain r > 2. We verify that globalExtSum( (1 , 0) ,S,N) 
agrees with (Ext^S^, iV)) >Q . 

15 : globalExtSum((l,0) ,S,N) 

05 = 

05 : S - module 

16 : prune truncate(0, Ext"l (truncate(2 , S~l) ,N) ) 

06 = 

06 : S - module 

Lastly, we show that (Extg(5>i, N)) >Q ^ 0. 

±7 : prune truncate(0, Ext"l (truncated , S"l) ,N) ) 

07 = cokernel {0} I OOOzOyOOOOxOOOwO I 

{OH z y i « I 
{0} I OzOOyOOOxOOOwOOO I 
{0} I OOzOOOOyOxOOOwOO I 

4 

o7 : S - module, quotient of S 

We conclude that the restrictions on r in theorem ^ are optimal. O 

Cohomology of Locally Free Sh eaves. Recall that, for a locally free sheaf C 
on X, there is an isomorphism | RH . propostion III. 6. 7]: 

Ext m (X ; £ v , Ox) = Ext m (X; O x ,C), 

where £ v = TComo x (£, Ox)- It follows from the equation (EH) that both of these 
global extension modules are isomorphic to H m (X, C). In particular, this gives two 
different methods for computing the cohomology of a locally free sheaf. In the next 
example, we compute the cohomology of a specific locally free sheaf by these two 
different approaches and compare the CPU time used. 
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Example 4.2. Consider the Veronese surface X embedded in P 5 . Following exam- 
ple 2.6 in Harris [JH], X is defined by the 2x2 minors of a generic 3x3 symmetric 
matrix. In Macaulay2, we construct the homogeneous coordinate ring R of X as 
follows: 

11 : kk = ZZ/32003; 

12 : S = kktUjVjWjXjjZ] ; 

13 : I = minors(2 ,genericSymmetricMatrix(S ,u,3) ) 

2 2 
o3 = ideal (-v +u*x, -v*w+u*y, -w*x+v*y, -v*w+u*y, -w +u*z, 

2 

-w*y+v*z, -w*x+v*y, -w*y+v*z, -y +x*z) 

o3 : Ideal of S 

14 : R = S/Ij 



Becau se X is smooth JH, exercise 14.13], the cotangent bundle fl is a loca lly free 
sheaf [R.H, theorem II. 8. 15]. Moreover, Eisenbud, Grayson and Stillman [ E-G-S| 
provide an algorithm called cotangentBundle which returns the module represent- 
ing the cotangent bundle. Using this, we create fl and f2 v : 

15 : Omega = cotangentBundle R; 

16 : OmegaDual = dual Omega; 

We establish the identity 0„> o Ext 1 ^; ft v , O x {y)) 0„> o Ext 1 (X; O x , Q(v )) . 
Evaluating the left hand side, we have: 

±7 : LHS = globalExtSum((l,0) .OmegaDual, R) 
o7 = cokernel {0} I z y x w v u I 

1 

07 : R - module, quotient of R 

The Macaulay2 output for the right hand side takes the form: 

18 : RHS = globalExtSum( (1 , 0) ,R, Omega) 

08 = cokernel {0} I z y x w v u I 

1 

08 : R - module, quotient of R 

Next, we use the benchmark function to produce an accurate timing for these 
techniques. The output is the number of seconds used by the CPU. 

19 : benchmark("globalExtSum((l,0) .OmegaDual, R)") 

09 = 0.51 

09 : RR 

ilO : benchmark("globalExtSum((1.0) ,R, Omega)") 

010 = 516.18 
olO : RR 

One immediately sees the significant difference in the length of two calculations. 
For this particular example, we can explain the difference by examining where the 
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first module is truncated. For the left hand side, no truncation is required since 
pd s (sR) =3 < 4 = n — I. Hence, evaluating globalExtSum( (1 , 0) , OmegaDual , R) 
involves finding the minimal free resolution of the i?-module Q and the S'-module 
sR. However, for the right hand side, we truncate at r = 2; thus calculating 
globalExtSum( (1 , 0) ,R, Omega) requires the minimal free resolution of i?-module 
i?>2 and the ^-module gfi. Now, the number of minimal generators of i?>2 is 
significantly larger than the number for $R and this leads to the difference in the 
amount of time used. From this example, we conclude that, at least for locally free 



sheaves, algorithm 3.5 and algorithm 3.6 are not the most efficient methods. O 



Duality Theory. For the third application, we illustrate Serre-Grothendieck dual- 
ity. Recall that if X is Cohen-Macaulay, closed subscheme of P™ of pure dimension 
d, then there exists isomorphisms Ext d_J (X; Q, lux) — H^(X,Q) V for < j < d, 
where Q is a coherent sheaf on X, lux is the dualizing sheaf and v denotes the dual 



vector space [RH, theorem III. 7. 6]. 



Example 4.3. Let X be the Del Pezzo surface o f deg ree 4 in P 4 ; X is a complete 



intersection of two quadratic hypersurfaces in P 4 [RH, exercise III. 4. 13]. We define 
the homogeneous coordinate ring R of X as follows: 

11 : kk = ZZ/32003; 

12 : S = kkCvjWjX.yjZ] ; 

13 : I = ideal (w*x, y *z) ; 
o3 : Ideal of S 

14 : R = S/I; 



By definition, we know that the dualizing sheaf lux is isomorphic to Ox{ — 1) [ RH 
remark III. 4. 7]. For simplicity, we let Q be the sheaf corresponding to the cokcrncl 
of a generic symmetric matrix. In Macaulay2 this appears as: 

15 : omega = R~{-1}; 

16 : G = coker genericSymmetricMatrix(R,v,2) 

06 = cokernel {0} I v w I 
{0} I w x I 

2 

06 : R - module, quotient of R 

Finally, we evaluate both sides of the duality equation. For j — 0, we have 

17 : EE2 = globalExt(2,G,omega) 

2 

07 = kk 

07 : kk - module, free 

18 : HH0 = Hom(sheaf Cohomology (0 , G) ,kk) 

2 

08 = kk 



08 : kk 



- module, free 
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In the case j = 1, we obtain 

i9 : EE1 = globalExt(l,G, omega) 
2 

o9 = kk 

09 : kk - module, free 

110 : HH1 = Horn (sheaf Cohomology (1 ,G) ,kk) 

2 

010 = kk 

010 : kk - module, free 

and, when j = 0, it is 

111 : EEO = globalExt(0,G, omega) 

011 = 

011 : kk - module 

112 : HH2 = Hom(sheaf Cohomology (2, G) ,kk) 

012 = 

This completes the example. O 



Global Extensions. As a final example, we compute a nontrivial global extension. 
We reminder the reader that a global extension of M. by Af is an exact sequence of 
Ox -modules: — ► J\f — ► £ — ► M — ► 0; two extension are equivalent if there is 
a commutative diagram 



0- 


->M — ► £ - 




* 










0- 


->A/" — >£' - 


-^M — 


-0 



The global extension modules derive their name from the following: the equivalence 
classes of global extensions are in bijective correspondence with Ext 1 (A; M., M) 
|G-H , section 5.3]. 



Applying theorem [l], there is a bijection between equivalence classes of global 
extensions and (Ext)j(M> r , JV)) Q for r sufficiently large. Moreover, the functor 
which takes a module to its associated sheaf is exact, so each global extension 

— >N — >£ — >M — >0 

corresponds to an exact sequence of i?-modules 

— > N — > E — ► M> r — > 0, 

again for r ^> 0. 



Example 4.4. Let X be the plane elliptic curve defined by x 3 + y 3 — z 3 . We will 
construct the unique nontrivial extension of Ox by Ox ■ We divide this computation 
into seven parts. 

(i) We create the homogeneous coordinate ring R of X. 
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11 : kk = ZZ/32003; 

12 : S = kk[x,y,z] ; 

13 : I = ideal (x~3+y~3-z~3) ; 

03 : Ideal of S 

14 : R = S/I 

04 = R 

04 : QuotientRing 

(ii) We calculate the k- vector space Ext 1 (X, Ox, Ox)' 

15 : globalExt(l,R,R) 

1 

05 = kk 

o5 : kk - module, free 

As Ext 1 (X; Ox, Ox) is one dimensional, there exists only one nontrivial 
global extension up to equivalence. 



(iii) We determine where to truncate the first module. Following algorithm 3.3 
we first calculate n and I. 

16 : n = (numgens S)-l, 1 = min(dim(R) , 1) 
06 = (2, 1) 

06 : Sequence 

We see that one must choose r > max {ai(sN) — i:l<i< pd s (siV)}. After 
computing F,(sN), we simply set r equal to the maximum: 

17 : sN = S'l/I 

07 = cokernel {0} I x3+y3-z3 I 

1 

07 : S - module, quotient of S 

18 : FF = res sN 

1 1 

08 = S < — S 

1 

08 : ChainComplex 

19 : r = (max degrees FF_1)#0-1 

09 = 2 

We introduce the modules M' and N to help avoid confusing notation. 

ilO : M' = prune truncate(2,R"l) 



olO = cokernel {2} lOyO-zOOOOx I 
{2} I y -z x I 
{2} I -z y x I 
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{2} I y x -z I 
{2} I -z x -y -z I 
{2} I x -y -z I 

6 

010 : R - module, quotient of R 

111 : N = R"l 

1 

011 = R 

oil : R - module, free 

(iv) We fix an exact sequence — ► K — ► P — ► M' — > 0, where P a free 
i?-module. This sequence is constructed from a free presentation of the R- 
modulc M', that is F 1 -A F — ► M' — > 0. 

112 : mu = presentation M' ; P = target mu; beta = map(M', P, id_P) ; 
il4 : K = image mu; alpha = map(P,K,mu); 

We verify that we have an short exact sequence. 

117 : {ker alpha == 0, image alpha == ker beta, image beta == M'} 
ol7 = {true, true, true} 

ol7 : List 

Now, this short exact sequence yields an exact sequence 

Hom fl (P, N) — > Rom R (K, N) — > Ext^(M', N) — > 0. 

Thus, an element in Ext fl (M, N) gives rise to an element 8 in Hornjj(if, N). 

(v) We choose a homomorphism 8: K — > N directly. 

118 : morphisms = Hom(K, N) 



image {-3} 


1 y 





-z 











x2 





{-3} 


1 y -z 














x2 





xz 


{-3} 


1 


y 





-z 











x2 


{-3} 


1 -z 





y 





X 











{-3} 


1 








X 


-y 








-z2 


{-3} 


1 





X 


-y 





-y z 


-z2 





{-3} 


1 


X 








-z 








-y2 


{-3} 


1 x 








-z 





-z2 


-y2 





{-3} 


1 x 


-z 











-y2 


-yz 





R - module, 


submodule 


of 


9 

R 











119 : choice = matrix{{l_R} , {0} , {0}, {0} ,{0} , {0} ,{0} , {0} , {0}} ; 

120 : theta = homomorphism map (morphisms , R~l, choice) 
o20 = lOyO-z 0000x1 

o20 : Matrix 

The module E is then the pushout of 8 and a; more explicitly the cokernel of 
the map tp : K — > P (B N where ^(a;) = (a(x), —9(x)) . 
(vi) We construct the module E. 
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i21 : D = P ++ N; 



i22 : psi = map(D, K, alpha I I -theta) 



o22 = 











-z 



y -z 
-z y 
y 
-z x 

-y 













x -z 



-z 



x -y 
-y -z 
-x 



o22 : Matrix 



123 : E = coker psi; 

(vii) We check that we have the short exact sequence — ► N — '"—t E — > M' — > 0. 

124 : iota = map(E, N, map(P,N,0) I I id_N) ; 

125 : phi = map(M' , E, id_P I map(P ,N, 0) ) ; 

126 : {ker iota == 0, image iota == ker phi, image phi == M'} 
o26 = {true, true, true} 

o26 : List 

Therefore, the i?-module E represents the unique nontrivial global extension. For 
completeness, we calculate the global sections of the sheaf associated to E. 

±27 : sheaf Cohomology (0,E) 
1 

o27 = kk 

o27 : kk - module, free 

As the space of global sections has dimension one, it follows that E represents an 



indecomposable rank 2 vector bundle [RH, theorem III. 2. 15]. O 
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